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MAPPING THE STABILITY FIELD OF JUPITER TROJANS. H.F. Levison, U.S. Naval 
Observatory, Flagstaff, AZ 86002, E.M. Shoemaker and R.F. Wolfe, U.S. Geological Survey, Flagstaff, AZ 
86001 

Jupiter Trojans are a remnant of outer solar system planetesimals captured into stable or quasistable 
libration about the 1 : 1 resonance with the mean motion of Jupiter. The observed swarms of Trojans 
may provide insight into the original mass of condensed solids in the cone from which the Jovian planets 
accumulated, provided that the mechanisms of capture can be understood [lj. As a first step toward this 
understanding, we have Undertaken to map the stability field of Trojans in the coordinates proper eccentricity, 
c p> and libration amplitude, D. To accomplish this mapping, we integrated numerically the orbits of 110 
particles with e p in the range 0 to 0.8 and D in the range 0° to 140°. Orbits of the Sun, the four Jovian 
planets, and the massless particles were integrated as a full N-body system, in a barycentric frame using a 
fourth order symplectic scheme [2, 3], Initial positions and velocities for the planets were taken from [4]. 

Our definitions of e p and D for Trojans are an extension of those presented in |l). We utilized the 
analytic work of firdi [5] to obtain 

e 2 = [ecosw - ey cos (w, ± w b )] 2 + [esinw - ey sin (wy ± w a )] 2 , 
where e is the osculating eccentricity, u is the longitude of perihelion, 

<2 0 = sin- 1 (#+^/ 2 + ^|#^), *„ = cos-Hi + &i 2 + ^ 4 ), 

and the subscript j refers to elements of Jupiter. The parameter l « 8.58 x 10“ 3 Z?, a result not explicitly 
given by £rdi. 

In our study of highly eccentric orbits we found that the longitude of the libration point is not fixed 
at ±60° from the mean longitude of Jupiter but is a function of c p . This is illustrated by inserts in Fig. 1, 
that show the trajectories of three particles in the rotating frame. The trajectories illustrated correspond 
to small values of D\ e p = 0.03, 0.3, and 0.7 in frames A, B, and C respectively. Small circles in Fig. 1 
represent measured values of the angle between the libration point and the mean longitude of Jupiter, <f>^ , 
as a function of e p \ the solid curve is a fourth order polynomial fit to the observations, 

<p L = 1.047 + 1.56e p + 0.901^ + 3.39e 3 - 3.71eJ. 

We calculate D as follows: Let x = | sin ((<£ — <f>L + x/3)/2)|, then 

where ^ = rj — tyy, <f> is the difference in mean longitude between the Trojan and Jupiter, m is the mass of 
Jupiter in solar masses, and rf is the mean motion of the object. Solutions for the minimum and maximum 
values of <f> can be found [6] from x min = sin (a/3)/3A x / 2 and x max = sin (a/3 + 120°)/3A 1 / 2 , where 
A = — mrij /2E eq and sin a = (3A) 3 / 2 . The libration amplitude D is <f>max - 

We integrated the orbits for time intervals corresponding to 1500, 15 000, and 150 000 Jupiter orbital 
periods. Test particles were started with orbits in Jupiter’s orbital plane. In a preliminary series of 
integrations carried out to 15 000 orbital periods, where the initial inclination of the orbits was 10°, we 
found no significant difference in the limit of the stability field from that obtained with t = 0. 

The distribution of orbits that are stable over the different intervals of time that we investigated is 
illustrated in Fig. 2, where the limits of stability for time intervals of 17 800, 178 000, and 1 780000 years are 
contoured. Orbits with surprisingly high initial t p are stable for nearly 2 million years, but no orbits with D 
exceeding about 110° are stable for more than 17 800 years. With increasing time the stability field shrinks 
in both c p and D . The limit of the main field of observed Trojan asteroids shown in Fig. 2 may represent 
the approximate limit of stability for a time interval of 4.5 billion years. 

One recently discovered Trojan, 1989 BQ, lies well outside the main Trojan field (c p = 0.22, D = 17°). 
We suggest that the dynamical lifetime of 1989 BQ may be of order 10 s years and that it has been captured 
late in solar system history. Capture may have occurred by diffusion, successively, from a free orbit to a 
horseshoe orbit to tadpole orbits of increasing stability. If one Trojan bright enough to be discovered by the 
usual methods of asteroid search has been captured late in solar system time, most of the known Trojans 
may have been captured by a similar process at an early time when the flux of Jupiter- crossing bodies was 
many orders of magnitude higher than at present [lj. 
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Figure 1: Position of libration point, 4>l, as a function of proper eccentricity, e p . 



Figure 2: Map of stability of Jupiter TVojans in proper eccentricity (e p ) and libration amplitude (17) space. 
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